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Realistic quantum mechanical systems are al- 
ways exposed to an external environment. The 
presence of the environment often gives rise to a 
Markovian process in which the system loses in- 
formation to its surroundings. However, many 
quantum systems exhibit a pronounced non- 
Markovian behavior in which there is a flow of in- 
formation from the environment back to the sys- 
tem, signifying the presence of quantum memory 
effects [IH5]. The environment is usually com- 
posed of a large number of degrees of freedom 
which are difficult to control, but some sophis- 
ticated schemes for modifying the environment 
have been developed |6j. The physical realization 
and control of dynamical processes in open quan- 
tum systems plays a decisive role, for example, 
in recent proposals for the generation of entan- 
gled states [7— 9|, for schemes of dissipative quan- 
tum computation [10J, for the design of quan- 
tum memories [11 J and for the enhancement of 
the efficiency in quantum metrology |12j . Here 
we report an experiment which allows through 
selective preparation of the initial environmen- 
tal states to drive the open system from the 
Markovian to the non-Markovian regime, to con- 
trol the information flow between the system and 
the environment, and to determine the degree of 
non-Markovianity by direct measurements on the 
open system. 

The standard approach to the dynamics of open quan- 
tum systems employs the concept of a quantum Markov 
process which is given by a semigroup of completely posi- 
tive dynamical maps and a corresponding quantum mas- 
ter equation with a generator in Lindblad form [T3l [14] . 
Very recently, a toolbox for the engineering of such quan- 
tum Markov processes in a multi-qubit system of trapped 
ions has been realized experimentally [15] and technolog- 
ical developments have also allowed experimental stud- 
ies of quantum correlations in open systems p~6j [T7] . 
Within a microscopic approach quantum Markovian mas- 
ter equations are usually obtained by means of the 
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Born-Markov approximation which presupposes a weak 
system-environment coupling and several further, mostly 
rather drastic approximations. However, in many pro- 
cesses occurring in nature these approximations are not 
applicable, a situation which occurs, in particular, in the 
cases of strong system-environment couplings, structured 
and finite reservoirs, low temperatures, as well as in the 
presence of large initial system-environment correlations. 
In the case of substantial quantitative and qualitative de- 
viations from the dynamics of a quantum Markov process 
one often speaks of a non-Markovian process, implying 
that the dynamics is governed by significant memory ef- 
fects. Quite recently important steps towards the devel- 
opment of a general consistent theory of non-Markovian 
quantum dynamics have been made which try to rig- 
orously define the border between Markovian and non- 
Markovian quantum evolution and to quantify memory 
effects in the open system dynamics [T8H2T] . 

The measure for quantum non-Markovianity con- 
structed in [19] is based on the idea that memory ef- 
fects in the open system dynamics can be characterized 
in terms of the flow of information between the open 
system and its environment. It has been used recently, 
e.g., to describe this information flow in the energy trans- 
fer dynamics of photosynthetic complexes [2j [4] , and to 
characterize memory effects of the dynamics of qubits in 
spin baths [5 . Here, we present the results of an experi- 
ment which enables through a careful preparation of the 
initial system-environment states and quantum state to- 
mography not only to control and to monitor the transi- 
tion from Markovian to non-Markovian quantum dynam- 
ics, but also the direct determination of this measure for 
quantum non-Markovianity. 

Quantum memory effects are quantified by employing 
the trace distance D{p\,p2) = §tr|pi — p2\ between two 
quantum states p\ and p^. This quantity can be in- 
terpreted as a measure for the distinguishability of the 
two states [22-24 r Markovian processes tend to contin- 
uously reduce the distinguishability of physical states, 
which means that there is a flow of information from the 
open system to its environment. In view of this inter- 
pretation the characteristic feature of a non-Markovian 
quantum process is the increase of the distinguishability, 
i.e. a reversed flow of information from the environment 
back to the open system. Through this recycling of infor- 
mation the earlier states of the open system influence its 
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FIG. 1: The experimental setup. Here, the abbreviations of 
the components are: HWP - half wave plate, QWP - quarter 
wave plate, IF - interference filter, QP - quartz plate, PBS - 
polarizing beamsplitter, FP - Fabry-Perot cavity, and SPD - 
single photon detector. 



later states, which expresses the emergence of memory 
effect in the open system's dynamics pT9^ [20] . 

On the basis of this physical picture one can construct 
a general measure for the degree of non-Markovianity 
of a quantum process given by some dynamical map & t 
which maps the initial states p(0) of the open system to 
the corresponding states p(t) = $> t (p(0)) at time t. The 
full time evolution of the open system over some time 
interval from the initial time to the final time T is then 
given by a one-parameter family ( £ = { < £t|0<£<T} 
of dynamical maps. We define the rate of change of the 
trace distance by <r(£, pi,2(0)) = ^D(pi(t), p2(t)). The 
measure «A/"(3>) for the non-Markovianity of the process 
is then denned by 

J\f($) = max / dta(t,p li2 (0)). (1) 

Pl, 2(0) J a>0 

Here, the time-integration is extended over all subinter- 
vals of [0, T] in which the rate of change of the trace 
distance a is positive, and the maximum is taken over 
all pairs of initial states. The quantity ([!]) thus measures 
the maximal total increase of the distinguishability over 
the whole time-evolution, i.e., the maximal total amount 
of information which flows from the environment back to 
the open system. 

In our experiment the open quantum system is pro- 
vided by the polarization degree of freedom of photons 
coupled to the frequency degree of freedom represent- 
ing the environment. The experimental setup is shown 
in Fig. [IJ An ultraviolet Argon-Ion laser is used to 
pump two 0.3mm thick BBO crystals cut for type I down 
conversion process to generate arbitrary pure two-qubit 
states. A fused silica plate (0.04mm thick and coated 
with partial reflecting coating on each side, with about 
85% reflection probability at 702nm) is used as a FP 
cavity. The cavity is mounted on a rotator which can be 
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FIG. 2: The frequency spectrum of the initial state for various 
values of the tilting angle 0. 

tilted in the horizontal plane. A 4nm (full width at half 
maximum) interference filter is placed after the FP cavity 
to filter out at most two transmission peaks. The corre- 
sponding interference filter in the other arm is lOnm. The 
polarization and frequency degrees of freedom are cou- 
pled in a quartz plate in which different evolution times 
are realized by varying the thickness of the plate. A po- 
larizing beam splitter together with a half- wave plate and 
a quarter- wave plate is used as a photon state analyzer. 

The half wave plate HWP 2 and the tilted FP cavity are 
used to prepare the initial one-photon states 1^1,2(0)) = 
|^i,2) ® |x), where 

l^i,2) = ±{\H)±\V)), (2) 

with \H) and \V) denoting the horizontal and the ver- 
tical polarization state, respectively. The environmental 
state |x) = J dujf(uj)\uj) involves the amplitude f(uj) for 
the photon to be in a mode with frequency lj, which is 
normalized as J duj\f(uj)\ 2 = 1. The form of the fre- 
quency distribution \f(uj)\ 2 and thus the initial state of 
the environment can be controlled by the tilting angle 
of the FP cavity. Figure [2] shows how determines 
the structure of the frequency spectrum and, thus, the 
environmental initial state |%). By changing the initial 
state of the environment in the experiment we modify the 
open system dynamics in a way which allows us to ob- 
serve transitions between Markovian and non-Markovian 
quantum dynamics. 

In the first version of the experiment, HWP1 is fixed 
at zero degree to generate a two-photon state. Pho- 
ton 2 is directly detected in detector SPD in the end 
of arm 2 as a trigger of photon 1. Photon 1 first passes 



3 



(a)" 




• = 6.0" 

■ = 7.83" " 

♦ = 2.5" - 
a0=1.5" 



effective path difference(unit of Ao) 




effective path difference(unit of Ao) 



FIG. 3: The trace distance (a) and the concurrence between 
the system and the ancilla (b) as a function of the effective 
path difference for four different values of the tilting angle 
0. The solid lines represent the theoretical predictions for 
cr = 1.8 X 10 12 Hz and Aw=1.6x 10 13 Hz. The effective path 
difference is equal to AnL and Ao = 702nm. The experimen- 
tal error bars due to the counting statistics are smaller than 
the symbols. 
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HWP2, preparing it in the state \(pi) or \if2)- The sub- 
sequent interaction between polarization and mode de- 
grees of freedom in the quartz plate is described by a 
quantum dynamical map <$> t acting on the open system, 
where the interaction time t is connected to the vari- 
able length L of the quartz plate by means of t = L/c. 
Finally, a full state tomography is carried out in detec- 
tor SPD in the end of the arm 1 to determine the po- 
larization state pi^(t) = 3>t (|^i,2) (^1,2 1) of photon 1. 
This allows the direct experimental determination of the 
trace distance D(pi(t), p 2 (£)) between the two possible 
one-photon states after a certain interaction time t con- 
trolled by the length L of the quartz plate. Experimental 
results for four different values of the tilting angle of the 
FP cavity are shown in Fig. [3^a). We clearly observe a 
crossover from a monotonic to a non-monotonic behavior 
of the trace distance as a function of time, i.e., a transi- 
tion from a Markovian to a non-Markovian dynamics of 
the polarization degree of freedom of the photon. 

The experimental results admit a simple theoretical 
analysis which is based on the fact that the time evolu- 
tion in the quartz plate may be described by the unitary 
operator U (t) which is defined by 



U(t)\\)®\uj) 



_ e in x ujt 



|A)<g>|o;), 



where n\ represents the refraction index for light with 
polarization A = H, V . The presence of the quartz plate 
thus leads to the decoherence of superpositions of po- 
larization states, which is due to a nonzero difference 
An = riy — Tin in the refraction indices of horizontally 
and vertically polarized photons. The corresponding dy- 
namical map <& t takes the form: 



\H)(H\ ^\H)(H\, 

\V)(V\ ^\V)(V\, 

\H)(V\ ^n*(t)\H)(V\, 

\V)(H\ ^K{t)\V)(H\, 



where the decoherence function is given by the 

Fourier transform of the frequency distribution, 



«(t) = J <M/(w)| 
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FIG. 4: The change of the trace distance and of the con- 
currence as functions of the tilting angle 0. The transition 
from the non-Markovian to the Markovian regime occurs at 
c± 4.1°, and from the Markovian to the non-Markovian 
regime at c± 8.0°. The positive values in the blue re- 
gions give directly the non-Markovianity measure Af(<&) of 
the process. The negative values in the grey area correspond 
to = 0, i.e., to Markovian dynamics. The error bars are 

due to the uncertainty of the tilting angle and the counting 
statistics. 



With the help of these relations it is easy to show that the 
trace distance corresponding to the initial pair of states 
([2| is equal to the modulus of the decoherence function, 



D( Pl (t), P2 (t)) = \K(t)\. 



(3) 



In the experiment the transition from Markovian to non- 
Markovian dynamics is observed through variation of the 
tilting angle of the FP cavity. As illustrated in Fig. [2j 
all frequency distributions are very well approximated by 
a sum of two Gaussians centered at frequencies uu^ with 
amplitudes Ak and equal widths a, which yields 



\K(t)\ = 

where A\ 



> 2 (Ant) 2 
1 + A 



y/l + A 2 + 2Acos(Au; • Ant), (4) 



1^4, A 2 = yqr4 and Auj = uj 2 
tilting angle of the cavity is relatively small and thus the 
distance Au between the peaks remains approximately 
constant. Therefore, the only relevant parameter control- 
ling the transition is the relative amplitude A. Equation 
Q yields an excellent approximation of the experimental 
data (see the continuous lines in Fig. [3]). 

The experiment also enables a direct determination of 
the measure for non-Markovianity ([!]). First, we note 
that the initial pair ^ is already optimal in the sense 
that it yields a maximal increase of the trace distance. 
The theoretical explanation for this fact is presented in 
the Supplementary Information. Second, in our experi- 
ment the increase of the trace distance is restricted to a 
single time interval (see Fig. [3|. The non-Markovianity 
measure ([T]) is thus obtained by determining the differ- 
ence of the trace distance between the first local mini- 
mum and the subsequent maximum. Our experimental 
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results are shown in Fig. [4] Increasing the tilting angle 
of the cavity decreases the relative amplitude A between 
the peaks in the frequency spectrum and thereby reduces 
the non-Markovianity of the process until a transition 
to Markovian dynamics occurs. Further increasing the 
tilting angle amplifies the relative amplitude again and 
brings the dynamics back to the non- Markovian regime. 

In Ref. [21 an alternative measure for non- 
Markovianity has been proposed which is based on the 
idea that a Markovian dynamics leads to a monotonic 
decrease of the entanglement between the open system 
and an isomorphic ancilla system, while a non-Markovian 
dynamics induces a temporary increase of the entan- 
glement. One can show (see Supplementary Informa- 
tion) that for the present experiment this measure co- 
incides with if one uses the concurrence [25] [26] as 
an entanglement measure. This fact leads to an alter- 
native and independent method for the measurement 
of the non-Markovianity by means of our experimen- 
tal setup. Fixing HWP1 to 22.5 degree, we generate a 
maximally entangled two-photon state. Photon 1 then 
passes the quartz plate and the composite final state is 
analyzed through two-photon state tomography. Exper- 
imental results are shown in Fig. J3^b) and Fig. |4j clearly 
demonstrating the equivalence of both measures for non- 
Markovianity. 

Our experiment clearly reveals the measurability of re- 
cently proposed theoretical measures for quantum non- 



Markovianity which yield important information about 
the type of quantum noise and about environmental 
properties, even when the environment is a complex 
system involving an infinite number of degrees of free- 
dom and is not directly accessible through measurements. 
Moreover, we have introduced a method for the control 
of the information flow between the open system and its 
environment, which opens the possibility of efficiently ex- 
ploiting memory effects in future quantum technologies 
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Supplementary Information 

Finding the maximizing pair of initial states 

For an arbitrary initial pair pi,2(0) the trace distance 
evolution is given by D(pi(t), p2(t)) = V / a 2 + \^(t)b\ 2 , 
where a = /^(O) — p^ify denotes the difference of the 
initial populations, and b = Pi 2 (0) — /?2 2 (0) the difference 
of the initial coherences. It follows that any growth of 



the trace distance is maximal for a = and \b\ = 1, in 
which case we obtain the formula Therefore, the pair 
of states given by Eq. ([2| is an optimal initial pair. More 
generally, all initial pairs of states are optimal which cor- 
respond to pairs of antipodal points on the equator of the 
Bloch sphere that represents the two-state state system. 



Equivalence of the two measures 

The initial state is the pure, maximally entangled 
system- ancilla state /?sa(0) = |^sa)(^sa| with |^sa) = 
(\HH) + | W)) . The dynamical map acts locally on 
the system part leading to the state 

PsaW = (*t®J)psA(0) 

= ±(\HH)(HH\ + \VV)(VV\ 

+K(t)*\HH)(VV\+K(t)\VV)(HH\y 

This state is diagonal in a basis of maximally entangled 
states (Bell state basis), the maximal eigenvalue being 
Pmax = ^(1 + I^WD- Therefore, the concurrence of psa(^) 
is given by C(psaW) = 2p max - 1 = [24 . Thus, we 

have the relation D(pi(t), p2(t)) = C(psa(£)) 5 showing 
that both measures yield identical results. 



